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The quantum phase transitions induced by a magnetic field are theoretically studied in a frustrated 
two-leg spin ladder. Using the density-matrix renormalization-group method, we find some magnetic 
phase transitions and plateaux in two different cases of strong and weak rung couplings. With the 
strong rung coupling, the three magnetization plateaux are found at 1/3, 1/2, and 2/3 due to the 
frustration. Those can be understood in terms of a quasi-spinon reconstructed from the singlet and 
the triplets of spins on a rung. The plateau at 1/2 corresponds to the valence bond solid of the quasi- 
spinons, while the plateaux at 1/3 and 2/3 can be associated with the array of quasi-spinons such 
as soliton lattice. This is different from the usual Bose-Einstein-condensatlon picture of triplons. 
Our results will be useful to understand magnetization curves in BiCu2P06. 
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I. INTRODUCTION 

Correlated many-body systems have rich variety of 
quantum phenomena such as Mott insulator, high tem¬ 
perature superconductivity, Kondo effect, fractional 
quantum Hall effect and so Interacting boson sys¬ 

tems also have many interesting physics such as super¬ 
fluidity accompanied by the Bose-Einstein condensation 
(BEC). A quantum spin system is useful for studying 
such an interacting boson system due to a mapping be¬ 
tween interacting spins- and bosons-systems^. In fact, 
the BEC in quantum anti-ferromagnet TlCuCb has been 
studied experimentally^i^ and theoretically^. Moreover, 
many quantum magnets have been examined so far—. 

The density of bosons in quantum magnets is tuned by 
an applied magnetic field, which induces the BEC corre¬ 
sponding to the long-ranged magnetic order. It can be 
seen in a magnetization curve. For example, in the dimer¬ 
ized quantum anti-ferromagnet TlCuCb, which has a gap 
in the spin excitation, the magnetic field has a critical 
value, Hci, where the gap is destroyed and the magneti¬ 
zation becomes finite. Such a behavior is explained by the 
BEC of triplons^i^i^, which are triplet states on a dimer 
and are associated with hard core bosons. Near H > Hci, 
the density of triplons is still dilute and will make some 
kinds of bound states^iiS. Increasing the magnetic field, 
H, the magnetization will become constant above a sat¬ 
uration field, Hc 2 - In Hci < H < Hc 2 , strong repulsion 
between triplons on a lattice can induce magnetization 
plateaux, which corresponds to Mott insulator phase of 
bosonsiip— . 

The triplon BEC mentioned above is accepted as 
a good starting point of understanding magnetization 
plateaux. In this paper, we proposed a new concept 
for magnetization plateaux, which is different from the 
triplon picture. Such a new concept works on a frus¬ 
trated two-leg spin ladder (F-2LSL). Calculating mag¬ 
netization curves of the F-2LSL by the density-matrix 


renormalization-group (DMRG) method, we find three 
fractional magnetization plateaux due to frustration in 
the case of strong rung coupling. We interpret the ori¬ 
gin of the plateaux by introducing a quasi-spinon con¬ 
structed by a singlet and a triplet of a spin pair on a 
rung. In contrast to the triplon BEC picture, the mag¬ 
netization plateaux correspond to the valence bond solid 
and the solitonic lattice of the quasi-spinons. As a realis¬ 
tic material containing the F-2LSL, BiCu2P06 (BCPO) 
has been actively studiedi^^— . Casola et al. have re¬ 
ported that a field-induced phase can be fit by a solitonic 
excitation of spinons originating from triplons instead of 
the triplon BEC pictureii, indicating a reconstruction of 
the quasi-particles and field-induced second-order phase 
transition. Our new concept of quasi-spinons will give a 
hint for fully understanding the field-induced phases of 
BCPO as well as other frustrated spin chains that exhibit 
fractional magnetization plateaujsii^— . 


The rest of the paper is organized as follows. In Sec.llll 
the Hamiltonian of the frustrated two-leg spin ladder will 
be defined, and a way to calculate the magnetization 
curves using the DMRG method will be explained. In 
Sec. mi numerical results of the magnetization curves, 
which depend on both the rung couplings and the frus¬ 
tration, will be shown in detail. We will formulate quasi¬ 
spin operators to explicitly derive an effective Hamilto¬ 
nian with the strong frustration and the strong rung cou¬ 
plings in an applied magnetic field. How to interpret the 
magnetization plateaux in terms of the effective Hamil¬ 
tonian is also discussed together with schematic pictures. 
Relation between weak and strong rung-coupling limits 
will be summarized in a table. Finally, summary and 
discussions will be given in Sec. IIVI 
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II. MODEL AND METHOD 


A model Hamiltonian of the F-2LSL with an applied 


magnetic field H is given by 


n = n\i+n± + nz 

( 1 ) 

with 


'H\\ = ’^v'Y^ 

77 = 1,2 j i—u,\ 

( 2 ) 

T-L_l = Jx Y^ ■ Sj^i, 

j 

(3) 


(4) 


j i=u,l 


where is the S' = 1/2 spin operator on the j site 

in the upper (lower) chain and its 2 component is 
There are three types of anti-ferromagnetic Heisenberg 
interactions, a nearest-neighbor coupling Ji and a next- 
nearest-neighbor coupling J 2 in the leg direction, and a 
nearest-neighbor coupling on a rung bond J±. In the 
limit of weak rung coupling, J± <C Ji, this model de¬ 
scribes decoupled two frustrated spin chains, while a non- 
frustrated spin ladder is obtained in the limit of weak 
frustration, J 2 <C Ji- Consequently, this model bridges 
between a frustrated spin chain and a non-frustrated spin 
ladder through Jj_ and J^- 

The preceding studies on this model have shown that 
there are two different phases as the ground state, i.e., 
the columnar-dimer and rung-singlet phases^^. In the 
columnar-dimer phase, the ground state is composed of 
two degenerated states with spontaneously-broken trans¬ 
lational symmetry 2 ^, while no degeneracy exists in the 
ground state of the rung-singlet phase. In the former 
case, spinon as a fermionic quasi-particle helps us to un¬ 
derstand the magnetic behavior, although the hard-core 
bosonic particle, triplon, gives a better explanation in the 
latter one. We note that this is one example of recon¬ 
struction of quasi-particles caused by the second-order 
phase transition with spontaneously-broken symmetry. 

We calculate the magnetization curves in both the 
weak- and strong-rung-coupling limits by using the 
DMRG method with an open boundary condition^^. 
First, we calculate the minimum energies Em with fixed 
magnetizations m = 0, l,...,Msat without a magnetic 
field. Then, the energy difference between Em and Em+i 
determines a magnitude of a magnetic field, 

Em-.m+l — Em+1 Em^ (5) 

at which the magnetization of the ground state changes 
from m to m + 1. Therefore, using Hm-.m+i and magnetic 
field the magnetiztion curve M(H) is caluclated by, 

-^sat 1 

M{H)= Y. m0{H-Hm-l:m)0{Hm-.ra+l-H), {6) 

m—1 

where the 0{x) is the Heaviside step function. 


III. NUMERICAL RESULTS AND QUASI-SPIN 
TRANSFORMATION 

Figures [ija)-(c) show that magnetization curves for 
rung couplings = 0 . 1 , 1 , and 10 , with a fixed frus¬ 

tration J2/J1 = 0.6 in a 72-rung system. In the weak 
rung-coupling limit, J±/Ji =0.1, as shown in Fig. [TJa), 
we can find a cusp singularity and a plateau at magne¬ 
tization ratio M/Msat = 1/3, which are also reported 
in a frustrated chaini^. With an applied magnetic field, 
a three-fold degeneracy of triplet excitation breaks and 
one state of them goes down to the singlet ground state. 
When the lowest energy of the triplets reaches down to 
that of the singlet ground state at El = Hd, the gap-to- 
gapless transition occurs and the magnetization becomes 
finite. A cusp singularity appears as a result of Lifshitz 
transition of Jordan-Wigner fermion^!. The strong frus¬ 
tration induces a one-third plateau originated from a re¬ 
constructed unit cell with the size three times larger than 
that of an original unit cell. Freezing spins with a mag¬ 
netic moment make a long-ranged order in the leg direc¬ 
tion together with deconfinement of spins between two 
legs (see Fig. [2ja)). Increasing the rung coupling up to 
J±/Ji = 1 in Fig. [IKb), the one-third plateau is weak¬ 
ened or disappears, although other plateaux emerge at 
one-half magnetization M/Mgat = 1/2 and two-third one 
M/Msat = 2/3. Moreover, with a strong rung coupling 
in Fig.jljc), we can find again the one-third plateau with 
some cusp singularities. The three plateaux at magneti¬ 
zation ratios M/Mgat = 1/3, 1/2 and 2/3 does not satisfy 
the condition proposed by Oshikawa et. ah^, which is 
given by, Q{S — m^) G Z, with spin S, periodicity of lat¬ 
tice Q {Q = 2 for our Hamiltonian), magnetization = 
MS/Mgat at plateaux, and the set of integers Z. In ac¬ 
cordance with the Oshikawa-Yamanaka-Affleck criteria, 
these are caused by spontaneously-broken symmetries of 
the plateau states in the present system. Since such a 
plateau state requires a confirmation with other meth¬ 
ods like the numerical calculation and a bosonization ap¬ 
proach one by one, these could be non-trivial, although 
potential plateaux and a possible case of spontaneously- 
broken symmetries have been discussed in several quan¬ 
tum spin models with similar characteristic a^^d^i^^i^^ . 

The frustration dependence of the magnetization curve 
with the strong rung coupling J±/Ji = 10 is also cal¬ 
culated with J2/J1 = 0.1, 0.4, and 0.7 in Figs. [Hd), 
(e), and (f), respectively. It is clear that cusp singu¬ 
larities and plateaux are caused by the frustration con¬ 
trolled by J2/Ji, even though the number of singulari¬ 
ties and plateau is different from the weak rung-coupling 
cases shown in Figs.[TKa)-(c). The spin freezing and long- 
ranged order with the one-third plateau state appear at 
J2/Ji = ■ In contrast to one-third plateau for weak- 

rung coupling, the periodicity of the freezing spins must 
be identical between the upper and lower chains because 
of the strong-rung coupling (see Fig. [2Kb)). 

In order to clarify the origin of the magnetization 
plateaux, a quasi-spin transformation^ is applied into the 
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Fixed to Ji! J\ = 0.6, 72 rungs 



Fixed to Ji / Jx = 0.1, 72 rungs 



FIG. 1. (Color online) Magnetization curves in a 72-rung F-2LSL. (a)-(c) Rung-coupling (J_i/Ji) dependence with a fixed 
frustration J 2 /J 1 = 0.6. (d)-(f) Frustration (J 2 /J 1 ) dependence with a fixed rung coupling Ji/J± = 0.1. The one-third plateau 
denoted by “1/3” and the cusp singularity denoted by “Cusp” in (a) correspond to those of a frustrated spin chain. Other 
plateaux denoted by “1/2” and “1/3” appear at Jx/Ti = 1.0 (b) and Ji/Jx = 0.1 (c). The 1/3 plateau reappears in (c) 
together with some cusp singularities. In (d) there are no plateau and cusp, but in (e) two cusp singularities and 1/2 plateau 
appear. Strong frustration at J 2 /J 1 = 0.7 causes three plateaux in (f). 


(a) weak rung-coupling limit (b) strong rung-coupling limit 



FIG. 2. (Color online) The schematic spin configurations for the plateau phases, (a) Weak rung-coupling limit and (b) strong 
rung-coupling limit, (c) Schematic configurations used in (a) and (b). In (a), up spins locally freeze with long-ranged order and 
total in each three neighboring spins in the upper (lower) chain equals to 1/2. In (b), similar quasi-up-spin freezing with 
long-ranged order occurs in the 1/3 and 2/3 plateau phases, although quasi-spins play a physical role instead of the real spins. 
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original Hamiltonian ([T} with a strong rung coupling. In 
this limit, the hard-core-bosonic picture is still applicable 
up to i? ~ Hci- However, the picture becomes worse in 
such a large magnetic field, because the symmetry of the 
triplets is no longer alive. Instead of the symmetry of the 
triplets, a new SU(2) symmetry constructed by a singlet 
and a triplet appears with a quasi-spin transformation 
given by. 


tL = ± * 




1 r 


rpZ _ 

“ 2 




(7) 


and 


= — 

J.m ^2 L 

. 1 

rjiZ _ 

^1“ “ 2 L 




-{Sl^ + Sl{) + S^^S-, + S-^S^, . ( 8 ) 


With the transformation, T operators obey the spin 
SU(2) algebra: {T+^,T-^} = 1 and = 0 for 

i ^ j and a = p, m. Therefore, the transformation makes 
us possible to consider another SU(2) quasi-spin model 
of two different energy scale discussed below. The origi¬ 
nal Hamiltonian with no perturbations (Ji = J 2 = 0) is 
rewritten as, 


1~L±^ “t“ '^z — ^ ^ ^ ^ ^ “t“ const., (9) 

j j 

where nj^a{= T^a + 1/2) is number operators of the 
quasi-spins. The effective chemical potentials for the 
quasi-spins are given by = J± — H and /rym = 
J_L (1 — nj^p) + H. With perturbations Ji ^ J 2 <C J_l, we 
consider magnetic fields |J_l —iJ| ~ Ji ~ J 2 <C J_l(~ H). 
Since the number operator rij^p is zero or one, the chem¬ 
ical potential ^ H is much greater than ^j^p ~ 0. 
Thus, the energy scales of quasi-spins split off by a strong 
magnetic field, where the low-energy physics is described 
by the quasi-spin of ‘p’. Since the magnetization process 
is obtained by the expectation value (n.j,p — as a 

function of H and 0, we can focus the discussion 

on the low-energy physics described by the quasi-spin 
of ‘p’. To neglect the high-energy physics, we consider 
the effective Hamiltonian in a Hilbert space projected by 
'Hes = VHV with V = Ilj/l^^i.m)- Here, the projected 
space has the expectation value of the number operator 
of the ‘m’ quasi-spin = 0. In this approximation, 

the up and down states of the quasi-spin of ‘p’ correspond 
a triplet and a singlet states on a bond, respectively (see 
also Fig. 2). In the projected Hilbert space, we finally 
find a frustrated spin chain again, where spins are sub¬ 
stituted with the quasi-spins as follows^ i^^i^^ , 


( 10 ) 


with 


eta—1,2 j 


jZfrpZ rpZ 


( 11 ) 


jxf 

I T~ -I-T“ T+ 1 

^ 2 ^ I’P i+'^’P J,p J+v,p> 


^'z = -H'Y.np^ 


( 12 ) 


where the XY- and Ising-components of quasi-exchange 
interactions are denoted by J^' = Jr, and J^' = Jj,/2 
with rj = 1,2, respectively. See the Appendix for more 
details, by which one can find possible extensions of the 
theory. This model has several differences from the orig¬ 
inal spin model in the weak rung coupling limit (see Ta¬ 
ble H]). One is anisotropy of the anti-ferromagnetic spin 
interactions, where every quasi-exchange interaction is 
XY-like, = 2J*'. In addition, we should be careful 
on the meaning of the magnetization of up and down 
quasi-spins, because the transformation maps the sin¬ 
glet state to down and the triplet to up quasi-spins. 
Thus, the total magnetization with g-filling of up spins 
is not 2{q — 1/2) but q. It must be also noted that the 
quasi-spins couple to the quasi-magnetic field given by, 
H' = H — J± — (Ji -I- J2)/2. This causes the zero ex¬ 
pectation value of number of up spins without the real 
magnetic field H = 0, that is, the singlet states pack in 
rung bonds in the ground state. 


With the help of the quasi-spin model, we can ap¬ 
proach the natures of the magnetization process in the 
strong rung-coupling limit. Firstly, the magnetization 
M/Msat = 1/2 with the strong rung coupling corre¬ 
sponds to M/Msat = 0 in the weak rung-coupling limit. 
The gapped ground state in the weak rung-coupling 
limit is understood by the Majumdar-Gosh state with 
spontaneously-broken translational symmetry, where the 
spin singlet dimers pack over the legs. In the same man¬ 
ner, the quasi-spin singlet dimers pack over the legs, al¬ 
though the magnetic moment M = -|-1 is set on a sin¬ 
glet dimer in Fig. [JKb). The 1/3 magnetization plateau 
with the weak rung coupling is equivalent to 2/3 (= 
l/3x 1/2-I-1/2) plateau in the strong rung-coupling limit. 
The freezing-spin structure and long-ranged order are 
also found in the strong rung-coupling limit, where the 
real spins of the upper and lower chains coherently freeze. 
Additionally, we can find a l/3(= —1/3 x 1/2 -|- 1/2) 
plateau and some cusps in the strong rung-coupling limit, 
which correspond to the M/Msat = —1/3 plateau and 
cusp singularities with low magnetic helds in the weak 
rung-coupling limit, respectively. 

Therefore, the plateaux in the strong rung coupling 
correspond to M/Mgat = ±1/3 and M/Mgat = 0 plateaux 
of the frustrated spin chain, although physical roles are 
played by the quasi-spins (see Fig. HKb)). Actually, we 


nes = v{n\^+n'z)v 
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TABLE I. Comparison between two different limits. 



Weak rung-coupling limit 

Strong rung-coupling limit 

Effective model 

Decoupled frustrated spin-4 chains Frustrated quasi-spin-4 chain 

Spin interaction 

Heisenberg 

XY-like 

Spin operator 


T 

M/Msat for g-filling of up spins 

2{q - 1/2) 

q 

Effective magnetic field 

77 

77' = 77 - Jx - (Ji -4 J2)/2 

Filling of up spins at 77 = 0 

1/2 

0 


confirm that the spins freeze coherently between the up¬ 
per and lower chains in the 1/3- and 2/3-plateau states, 
and that quasi-spin correlation rapidly decays in the 1/2- 
plateau state. Moreover, we can see good correspondence 
between Fig. [TJa) with the weak rung-coupling and the 
quarter area at the upper right in Fig. [TJc) with the 
strong rung-coupling including a cusp singularity. Since 
the filling of the quasi-up-spins with the strong rung cou¬ 
pling corresponds to that in the ground state without a 
magnetic field in the weak rung coupling limit, the mag¬ 
netization curve in Fig. [Ija) is folded in half in Fig.jljc). 
Thus, we conclude that the reconstruction of the quasi¬ 
particles results in the emergence of the isostructural 
magnetization curves in two different limits. 


IV. SUMMARY 

In summary, we have studied the magnetization curve 
of the frustrated spin-ladder system using the DMRG 
method, and the three magnetization plateaux are found 
at M/Msat = 1/3, 1/2, and 2/3 due to frustration. In 
the limit of strong-rung coupling, the magnetization near 
the field-induced phase transition [H ^ Hd) can be un¬ 
derstood by the triplons on rungs as hard-core bosons. 
On the other hand, the plateaux in Hd < H < Hc 2 
can be explained by the quasi-spinons, which are recon¬ 
structed from the singlet and the triplet of spin pairs 
on a rung. The plateau at M/Mgat = 1/2 corresponds 
to the valence bond solid of the quasi-spinons, while the 
plateaux at 1/3 and 2/3 can be associated with the array 
of quasi-spinons such as soliton lattice. This is different 
from the usual BEG picture of triplons. The magneti¬ 
zation curves around 1/3 and 2/3 look like to be folded 
in half with same characteristics. This is also naturally 
explained using the effective model of quasi-spinons with 
spontaneously-broken symmetry. 

Our results will be useful for frustrated quantum spin 
ladder systems such as BGPO. Concerning BCPO, it is 
noted that an anisotropic interaction, i.e., Dzyaloshinsky- 
Moriya (DM) interaction, is reported by the inelastic neu¬ 
tron scattering study^S. The DM interaction also change 
the usual BEG picture of triplons on the magnetization 
process in the spin ladder system^S. The transition will 
be modified such as the second order one due to the DM 
interactions, which lead to linear terms of triplons. When 
the DM interaction is smaller than the rung coupling, our 
picture of the quasi-spinons on the magnetization curve 


will not be changed. Effects of such an interaction on the 
magnetic plateaux will be clarified in the near future. 
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Appendix A: Derivation of effective Hamiltonian 
and qnasi-spin operators 

In this section, we present a derivation of effective 
Hamiltonian with a quasi-spin transformation in the 
strong rung-coupling limit, Ji/J_l <C 1 and J 2 /J 1 . “C 1. 
With a finite magnetization in this limit, a quasi-spin 
transformation and a reduced Hamiltonian are useful to 
understand the ground state and the low-energy physics. 
In order to obtain the reduced Hamiltonian, we start with 
two spin Hamiltonian on j-th rung as follows, 

= 

— -t- (rij^u + nj^\) + ^ H, (Al) 

with a Jordan-Wigner transformation of spin operators, 

= (A2) 

and the number operators Uj up) = = 

”*"5" As it is well-known, these Jordan-Wigner 
operators obey anti-commutation relation, {dj,i,dj f.} = 




6 


6i^k and = 0- To diagonal¬ 

ize this Hamiltonian, we can use a bonding and an 
anti-bonding operators for create and annihilate oper¬ 
ators of Jordan-Wigner fermions, dj^h = {dj,u + dj^i)/\/2^ 
dj,a. = (dj,u — dj^{)/'/2- With the number operator of the 
bonding and an anti-bonding Jordan-Wigner fermions 
rij b(a)j the rung Hamiltonian is rewritten as, 

= nj.aK'.b - 1) - + nj,b - 1) -f (A3) 

We can obtain a quasi-spin transformation Eqs. o and 
dSl) with an inverse Jordan-Wigner transformation as fol¬ 
lows, 


and 


t; 




T7 




rpZ 


1 


(A5) 

These operators also satisfy the spin SU(2) algebra for 
themselves, and commutate each other. The quasi-spin 
operators rewrite the rung Hamiltonian as. 


•^rung 






= P 





where effective quasi-magnetic fields. 





(A7) 


= id^ 


J^P 


'j,b’ 


Tjp — idj^b: 


^j\p — 2 ’ 


(A4) 


With the quasi-spin operators, the leg Hamiltonian be¬ 
tween j-th and (j + L)-th rungs is given by. 




z=u,l 

_ f'Y'Z _ rpZ \ (rpZ _ rpZ \ _ ('T'^ T''^ _ T‘~ T'~ \ i ^ 

2 \j^P j,m) j+L,p \-^j,p-^j,m j,p-^ j,mj j+L,p-^ j+L,m j+L,p-^ J 


^91 ^j^P^j+L,p 


COS 


L2 


(rpZ _ rjiz \ 


-L 7^+ T~ 
' j,Ta 


COS 


cos 




-T,.„,T,+i,pCos 


■(^;-p-^w 


H.c. ^ 


(A8) 


If we consider small quasi-magnetic field for Tj^-p spin, 
namely |iJp/T±| -C 1, the quasi-magnetic field for 
spin is much larger than |iJp|, namely |iT,,m| > ^ » 
|iJp|. To deal low-energy physics, we can project out 
the high-energy states, that is, quasi-up-spins of 
opeartors. With the projection operator given hy V = 
rij (Tym ~ effective Hamiltonian is obtained as 

'PTL'P = TLeS- Since the original Hamiltonian © is com¬ 
posed by a sum over the rung and leg Hamiltonians, 
and the quasi-spin transformation gives us 

the effective Hamiltonian Eq. m. 

We note that the quasi-spin operators can be written 
by singlet and triplets conhgurations on j-th rung as fol¬ 
lows, 

n, = - \n,{t-\, - \s)j{s\,). 

^i!p = a-nd H.c., (A9) 

and 

- \s),{s\,), 

Tym = \t~)j{s\ 3 and H.C., (AlO) 


where \t°‘)j {a = ±,0) and |s)j denote the triplets and 
singlet states on j-th rung, respectively. When the quasi- 
up-spin states of operators are projected out, the 
quasi-up-spin and quasi-down-spin of Tj^p operators ap¬ 
proximately correspond to the triplet state and the 
singlet state |s)j, respectively (see Fig.[3|). 



FIG. 3. (Color online) Schematic relationship between the 
quasi-spin operators and the cofigurations of the singlet and 
triplets on J-th rung. When the quasi-up-spin states of Tj^m 
operators are projected out, the left-down side is prohibited. 
Thus, the quasi-up-spin and quasi-down-spin of Tj.p operators 
approximately correspond to the triplet state \t^}j and the 
singlet state |s)j, respectively. 
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